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NUMERICAL  DECONVOLUTION  OF  PHOTOELECTRON  SPECTRA  OF  LIQUIDS 

James  B.  Flanagan* 

Department  of  Chemistry,  New  York  University 
4 Washington  Place,  New  York,  N.  Y.  10003 

Abstract 

Described  is  a method  for  the  deconvolution  or  removal  of  spectral 
broadening  due  to  electron  interactions  with  solvent  in  liquid-phase 
ultraviolet  photoelectron  spectroscopy.  Certain  differences  between  this 
usage  of  deconvolution  and  its  usual  applications  in  pes  are  noted. 
Application  of  the  technique  is  illustrated  using  synthetic  and  experimental 
data.  Possible  uses  include  recovering  the  "original"  photoemission  spectrum 
before  broadening  by  liquid  as  well  as  studying  the  mechanism  of  energy  loss 
in  solvent.  Effects  of  unscattered  electrons  on  the  interpretation  of  this 
procedure  are  discussed. 

INTRODUCTION 

Although  a theory^  of  the  spectral  broadening  due  to  interactions  of 

photoemitted  electrons  with  liquid  in  bulk  phase  has  been  available  for  some 

time,  and  its  validity  recently  confirmed  on  a more  mathematically  rigorous 
2 

basis,  there  has  remained  a gap  in  the  application  of  this  theory  to  the 
detailed  interpretation  of  liquid-phase  photoelectron  spectra.  Up  until  now, 
interpretation  of  these  electron  energy  distribution  curves  (EDC's)  has  been 
made  on  the  basis  of  semi -quantitative  measures  such  as  threshold  energy  for 
emission,  or  by  more  sophisticated  measures  such  as  superposition  of  EDC's  or 
second  derivative  curves  (SDC's)  which,  although  able  to  reveal  something 
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about  the  detailed  spectra,  rely  on  certain  simplifications  or  assumptions  for 
their  justification  (although  SDC  curves  may  be  corrected  empirically,  band  by 
band,  to  approximate  the  results  of  the  deconvolution  method  described  here). 
The  present  technique  Is  related  to  the  deconvolution  procedures  of  other 
spectroscopies,  and  can  be  used  to  remove  the  effects  of  spectral  broadening 
In  liquid,  at  least  Insofar  as  these  effects  are  adequately  described  by  the 
theory  of  energy  loss  In  bulk  liquid. 

THEORY 

Liquid-phase  pes  spectra  of  electrons  originating  In  the  "bulk"  of 
solution  can  be  described  by  a convolution  Integral, 

f®(T)  - J"f^x)B(T,x)dx  (1) 

where  f^’CT)  Is  the  measured  partial  current  of  electrons  of  kinetic  energy 
T,  f^  Is  the  "true"  electron  current  at  the  point  of  emission,  before 
scattering  by  liquid,  and  B(T,x)  Is  the  convolution  function.  Equation  (1)  Is 
a slightly  generalized  form  of  Eq.  (1)  of  Ref.  3,  which  has  been  modified  here 
to  take  into  account  the  dependence  of  the  convolution  function  on  the  kinetic 
energy  of  the  electron  upon  emission.  A more  usual  application  of 
deconvolution  techniques  is  the  removal  of  broadening  due  to  nonmonochromatic 
photon  sources  or  nonideal  electron  energy  analysis.  In  this  case,  B is 
usually  given  as  a function  of  (T  - x).  There  are  numerous  methods  of 
obtaining  approximate  solutions  to  Eq.  (1),  many  of  them  Involving  some  degree 
of  Implicit  smoothing.  It  was  found  to  be  expedient  to  solve  Eq.  (1)  by  a 
so-called  direct  method.  In  which  deconvolution  of  a given  f*"(T),  as  a 
digitally  tabulated  array,  yields  a unique  f^{T)  array,  according  to  the 
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where  and  T"  are  Interpreted  as  data  vectors  and  B the  matrix 

3 4 

corresponding  to  the  convolution  function.  * It  can  be  shown  that  the 
recursive  method  described  below  for  the  solution  of  Eq.  (1)  corresponds  to 
the  solution  given  In  Eq.  (2).  A major  drawback  to  the  direct  method  Is  the 
problem  of  noise  amplification.  Use  of  Fourier  Transform  smoothing 

e g 

techniques  * were  absolutely  necessary  for  the  experimental  data  before 
deconvolution.  More  details  of  the  Fourier  Transform  method  will  be  given 
below. 

In  the  theory  of  Nemec,^  the  Instantaneous  loss  of  energy  from  energetic 
electrons  to  solvent  molecules  In  the  course  of  random  walk  to  the  surface 
where  they  are  collected  and  energy  analyzed.  Is  described  In  a single 
equation  of  the  form, 

d(T^)/dt  - - av"  (3) 

where  T^^  Is  the  electron  energy  In  liquid,  t Is  time,  a Is  a constant,  v Is 
electron  velocity  (T^^  » 2 m^v  ),  n Is  a constant  which  depends  on  the 
physical  process  of  energy  loss  to  solvent. 

If  an  electron  becomes  thermal Ized  (T^^  0)  before  reaching  the  surface. 

It  is  not  counted  by  the  energy  analyzer.  Electrons  originating  from  deeper 
within  the  solution  are  less  likely  to  reach  the  surface  and  undergo  more 
energy  loss,  on  average,  than  do  electrons  originating  nearer  the  surface. 
Taking  these  factors  Into  account  In  the  random  walk  problem,  the  following 
equation  for  the  EDC  of  electrons  reaching  the  surface  can  be  derived.^ 
From  Eq.  (16)  of  Ref.  1,  generalized  for  any  value  of  n 3,  this  1 s: 

dN/dT^  - (QCF^/K2)(2x/3ff)^/2  J "D(Tg) 

' ^s 


(4) 
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where  dN(T^)/dTg  1s  the  EDC  function,  Q is  the  cross  section  for 
photoelectron  emission,  C is  emitter  concentration,  F is  photon  flux,  K is  a 
normalization  constant  (see  Ref.  1,  Eq.  (5)),  x is  mean*square  random  walk 
step  length,  D(Tg)  is  the  initial  distribution  of  electrons  as  a function  of 
energy  of  photoemission  (this  is  the  unknown  function  which  we  wish  to  find), 
Tg  is  kinetic  energy  of  electron  at  the  instant  of  photoemission,  is 
kinetic  energy  of  electron  at  the  surface  of  the  liquid. 

Combining  the  constant  terms  into  K^,  we  obtain 

^s 

Transposing  this  into  finite  difference  form: 

N 

E(T  ) = K,  I D(T  + jAT)((T  + jAT)P  - (6) 

s i s s s 

where  E is  the  discrete  EDC  data  set  tabulated  at  intervals  of  aT,  N is  taken 
sufficiently  large  that  0(1^  + jAT)  ■ 0 for  all  j ^ N,  p is  (3-n)/2. 

Given  that  E(T^)  is  known  for  all  T^,  it  is  possible  to  calculate 
D(Tg)  recursively: 

E(Tj)  » Ddj  + Ll)  K^((T5  > aT)P  - tP)^/2^T 
N 

+ K,  I D(JaT  + T )((T  + jAT)P  - tP)1'2aT  (7) 

1 J.2  S S s 


whence. 


Edj)  - I D(Tj  ♦ jAT)((T^  + jAT)P  - tP)^/^ 


D(Ts  + AT) 


Ki((T5  + aT)P  - 


(8) 


5 

Thus,  given  that  for  sufficiently  large  T^,  D(Tg)  Is  Identically  zero, 

it  is  possible  to  use  Eq.  (8)  recursively  to  calculate  the  array  of  D(Tg). 

Although  It  has  not  been  shown  explicitly  In  the  equations,  D(Tg)  Is  j’ 

tabulated  at  Intervals  of  AT,  from  large  T„  to  T„  * aT.  For  the  first 

9 9 

nonzero  value  of  E(Tp,  Eq.  (8)  becomes: 

D(T^  + AT)  . E{T^)/K^{{T^  + AT)P  - tP)AT  (9) 

The  constant  Is  automatically  taken  Into  this  value  of  D and  all 

subsequent  values  In  the  D array.  This  has  two  effects:  (1)  It  Is  not 

necessary  to  calculate  a priori,  which  would  be  nearly  Impossible,  and 
(2)  Kj^  will  differ  In  numerical  value  when  different  values  of  p (or  n)  are 
used  In  the  deconvolution,  making  the  results  appear  to  different  scale  (there 
Is  also  a difference  of  scale  between  deconvolution  and  SOC).  This  recursive 

.. 

method  corresponds  to  the  exact  solution  of  Eq.  (2). 

i 

RESULTS 

Simulated  EDC  data.  In  order  to  check  the  correspondence  of  the  | 

deconvolution,  Eq.  (8),  to  the  generating  function,  Eq.  (6),  synthetic  EDC  | 

curves  were  generated  from  a known  D function.  To  test  the  ability  of  the 

deconvolution  to  work  on  functions  of  arbitrary  complexity,  the  D function  was 
chosen  to  be  portions  of  two  Inverted  parabolas  superimposed  on  a flat  zero 
baseline.  Use  of  such  a function  enables  testing  of  the  ability  of  the 

deconvolution  to  separate  closely  spaced  and  rather  narrow  bands.  Figure  1, 
curve  A shows  the  D function  which  was  assumed  In  all  the  simulations  to 
follow.  Curve  B shows  the  EDC  calculated  from  A using  Eq.  (6)  assuming  p > 1 
(corresponding  to  dT^/dt  » - av).  Without  the  Introduction  of  random  noise 
or  Fourier  Transform  snoothing,  and  assuming  the  same  value  of  p used  In  the 

EDC  generation,  the  D function  was  reproduced  exactly. 
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Contrast  this  behavior  with  the  (uncorrected)  SDC  of  the  same  data  shown 
in  Fig.  2.  As  shown  previously  for  simulated  Gaussian  distribution,^  there 
is  a shifting  toward  lower  kinetic  energy,  broadening  of  the  peak,  and  a 
rising  baseline  toward  lower  energies.  Nevertheless,  the  SDC  does  seem  to  be 
capable  of  good  performance  with  respect  to  the  difference  in  energy  between 
successive  peaks,  and  reproduces  the  qualitative  aspects  of  the  spectrum 
rather  well. 

In  order  to  test  the  ability  of  the  EDC  deconvolution  technique  to 
discriminate  the  energy  loss  law,  Eq.  (3),  EDC  data  with  a certain  value  of  p 
were  deconvoluted  assuming  a different  value  of  p.  In  Fig.  3,  the  original 
data  were  generated  assuming  p » 1 and  fit  with  Eq.  (8)  assuming  p > 2.  In 
this  case,  there  is  a rising  portion  at  low  kinetic  energies  similar  to, 
though  not  as  severe  as  the  SDC  case.  The  relative  magnitude  of  the  two  peaks 
is  affected,  but  the  peak  positions  are  not.  Deconvolution  of  real  data 
assuming  a value  of  p which  is  too  large  would  therefore  be  a difficult 
situation  to  detect  because  the  systematically  high  baseline  could  not  be 
differentiated  from  real  emission. 


Figure  4 shows  the  situation  in  which  the  value  of  p assumed  in  the 
deconvolution  is  smaller  than  actually  used  in  the  generation  of  the  EDC. 
Here  we  see  that  areas  of  negative  emission  are  calculated.  Since  this  cannot 
happen  in  a real  physical  situation,  this  effect  might  be  thought  to  be 
diagnostic  of  a choice  of  p which  is  smaller  than  physical  reality.  Actually, 
disproportionate  contributions  from  unscattered  electrons  can  also  produce 
this  effect.  Again  we  see  that  the  relative  heights  of  the  two  peaks 
(initially  equal)  is  affected.  The  peak  energies  are  reassuringly  accurate. 

The  results  seen  in  Fig.  4 point  to  a crude  but  possibly  usable  method  of 
determining  the  value  of  p (i.e.,  the  energy  loss  law)  in  a real  physical 
system.  Necessary  would  be  an  emitter,  preferably  a pure  liquid,  which  is 
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(known  from  gas-phase  data  to  have  a region  of  very  low  (near  zero) 
photoelectron  emission  at  higher  binding  energies  than  the  first  emission 
band.  One  would  have  to  choose  such  a liquid  on  the  basis  of  a gas-phase 
spectrum  showing  an  emission  gap  at  least  2-3  eV  wide  (and,  of  course,  having 
I the  other  properties  of  a good  liquid-phase  pes  prospect).  By  taking  very 

i- 

j accurate  EDC  data  using  a source  of  low-energy  photons  to  minimize  the  yield 

\ 

of  unscattered  electrons,  the  value  of  p could  be  found  at  which  the  results 

!'■ 

of  deconvolution  yield  zero  emission  in  the  gap.  In  any  case,  a lower  limit 
on  p might  be  obtained  by  this  procedure.  Unfortunately,  as  shown  below,  the 
' existence  of  unscattered  electrons  may  make  such  a determination  of  the  energy 

loss  law  difficult  or  impossible. 


' 


Application  to  Experimental  Data.  In  order  to  test  the  applicability  of 
the  deconvolution  method  to  real  data,  EDC  data  for  ethylene  glycol  and 

Q 

N-methyl aniline  reported  previously  were  used.  The  data  were  read  from 
paper  tape,  where  they  had  initially  been  stored,  into  a PDP-11/34  laboratory 
computer  where  they  were  stored  on  floppy  disk.  The  data  were  originally 
taken  as  total  current  curves  at  a retarding  potential  grid  and  had  to  be 
differentiated  once  (with  respect  to  retarding  potential)  to  yield  the  EDC. 
Smoothing  at  this  stage  was  not  necessary.  The  zero  of  potential  was 

O 

determined  by  the  capacitance  technique.  These  capacity  curves  had  been 
recorded  on  paper  tape  at  the  same  time  as  che  emission  data.  All  EDC  and 
deconvolution  data  shown  are  already  corrected  for  this  offset  potential. 

Upon  going  from  the  EDC  to  either  the  SDC  or  deconvolute,  smoothing  was 
found  to  be  absolutely  necessary.  The  Fast  Fourier  Transform  technique^’® 
with  the  sigma  smoothing  factors  described  by  Lanczos^  was  used.  In  this 
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method,  the  EDC  spectrum  is  first  extended  with  baseline- level  currents  to 
obtain  cyclic  continuity  required  for  Fourier  transformation  and  to  yield  a 
data  set  numbering  a power  of  2,  ;n  the  cases  below,  512  points.  The 
transform  of  a real  array  of  512  data  points  yields  an  array  of  256  complex 
"frequency"  space  Fourier  coefficients.  Smoothing  is  accomplished  by 
truncating  this  array  of  coefficients  so  as  to  eliminate  higher  frequency 
components  in  which  noise  resides.  The  sigma  coefficient  method  is  one  way  of 
accomplishing  a "soft"  cutoff  which  tends  to  minimize  unwanted  oscillations 
near  the  cutoff  frequency  in  the  smoothed  data.  The  choice  of  cutoff  was  made 
by  visually  observing  the  tradeoff  between  noise  in  the  deconvolute  and  the 
loss  of  detail  due  to  excessive  smoothing.  In  the  case  of  the  deconvolution 
of  N-methyl aniline,  below,  special  care  was  needed  because  of  the  many  closely 
spaced  bands.  It  was  always  possible,  nevertheless,  to  choose  a range  of 

cutoffs  for  which  the  peak  energies  were  well  defined  and  unchanging,  and  the 

noise  was  not  excessive.  In  the  cases  below,  the  smoothing  was  effected  by 
truncating  the  components  above  about  the  40th  of  the  256  complex 
coefficients.  This  is  roughly  equivalent  to  a moving  smoothing  function  with 
a window  0.2-0. 5 eV  wide. 

Ethylene  Glycol.  Figure  5 shows  a portion  of  the  EDC  of  liquid  ethylene 
glycol  as  a function  of  kinetic  energy,  corrected  for  the  surface  potential  as 

O 

described  above.  The  data  were  taken  as  described  with  a He(I)  resonance 

lamp  photon  source  of  21.2  eV.  "Binding  energy"  can  be  calculated  as  (21.2  - 

T)  eV.  For  ethylene  glycol,  the  EDC  is  monotonically  increasing  (from  right 
to  left).  Thus,  it  should  be  possible  to  find  a value  of  the  power  law  for 
which  negative  emission  is  not  observed.  It  turns  out  that  the  choice  of  a 
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power  loss  law  coefficient  sufficiently  large  to  eliminate  the  areas  of 
"negative"  emission  seen  in  the  deconvolute.  Fig.  6 for  p * 1,  is  impossibly 
large,  p » 5 or  6,  with  severe  distortion  of  the  resulting  deconvolute.  The 

Q 

conclusion,  reinforcing  the  earlier  observations  of  Nemec  ^ ^.,  is  that 
there  exists  a strong  component  cf  unscattered  electrons  in  the  EDC  spectrum. 
The  deconvolution  method  appears  to  be  superior  to  the  method  of  empirical 
correction  of  SDC's  in  that  it  is  able  to  demonstrate  (by  the  existence  of 
negative  areas  in  the  deconvolute)  the  presence  of  un scattered  electrons  even 
when  the  EX  is  monotcoic. 

N-Methyl aniline.  This  material  exhibits  more  obvious  effects  of  emission 
of  unscattered  electrons  than  does  ethylene  glycol.  Figure  7 shows  the  EDC  of 
data  recorded  previously®  at  -44°C,  using  a He(I)  photon  source.  The 
presence  of  unscattered  electrons  is  indicated  by  the  maxima  in  the  high 
kinetic  energy  peaks.  Deconvolution  of  such  an  EDC  (Fig.  8)  indicates  even 
more  pronounced  negative  peaks  in  the  result  than  in  the  case  of  ethylene 
glycol. 


! 


i 


t 

i 

1 


[ 1 

1 • 


CONCLUSION  I 

It  has  been  shown  to  be  possible  to  perform  deconvolution  directly  on  ■ ■ 

liquid-phase  EDC  data  without  relying  on  corrected  SDC  analysis  assuming  j 

Gaussian  D function.  The  problem  of  the  overabundance  of  unscattered  i 

electrons  is  seen  to  be  a major  drawback,  however.  The  attempt  at  correction  j 

for  unscattered  electrons  by  superposition  of  corrected  SDC  derived  Gaussian  | , 

f « 
. < 

emission  from  bulk  with  a variable  band  due  to  unscattered  electrons  was  shown  | 

previously  to  be  effective  in  fitting  the  data,®  though  the  physical  i 
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meaningful ness  of  this  procedure  may  be  called  Into  question.  An  empirical 
modification  of  the  bulk  emission  theory  which  gives  more  weight  (as  a 
function  of  energy)  to  electrons  originating  near  the  surface  might  provide  a 
solution  to  the  problem,  L.^ugh  without  much  additional  physical  insight  and 
requiring  the  Introduction  of  extra  adjustable  parameters.  Ultimately,  what 
would  be  desired  would  be  a unified  theory  of  photoelectron  emission  from 
liquids,  with  detailed  consideration  of  the  discrete  nature  of  the  liquid 
surface  layer. 
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CAPTIONS  TO  FIGURES 


Figure  1.  Curve  A - Initial  distribution  D(Tg);  Curve  B - EDC  calculated  by 
Eq.  (6),  as  function  of  kinetic  energy.  Both  functions  were  calculated  as 
discrete  functions,  the  actual  tabulated  values  shown  as  dots  on  the  curves. 

Figure  2.  Results  of  taking  SDC  (derivative)  of  Curve  B in  Fig.  1.  Discrete 
points  not  shown. 

Figure  3.  Results  of  deconvolution  assuming  p « 1 for  calculation  of  EDC  (Eq. 
(6))  and  deconvoluted  assuming  p » 2 (Eq.  (8)). 

Figure  4.  Results  of  deconvolution  assuming  p * 2 for  calculation  of  EDC  and 
p » 1 for  deconvolution. 

Figure  5.  Portion  of  EDC  for  ethylene  glycol  after  Fourier  Transform 
smoothing.  Data  from  Ref.  8. 

Figure  6.  Deconvolution  of  glycol  EDC  data  of  Fig.  5.  p assumed  « 1.  Little 
difference  in  results  was  observed  for  p * 0 or  2. 

O 

Figure  7.  Smoothed  EDC  of  N-methyl aniline. 

Figure  8.  Deconvolute  of  EDC  in  Fig.  7.  Variation  of  p fro^ii  value  of  1 used 
here  had  little  effect  on  extensive  areas  of  "negative  emission." 
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